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ON THE INVERSE PROBLEM OF FINDING COSMIC 
STRINGS AND OTHER TOPOLOGICAL DEFECTS 

MATTI LASSAS, LAURI OKSANEN, PLAMEN STEFANOV, 

AND GUNTHER UHLMANN 


Abstract. We consider how microlocal methods developed for 
tomographic problems can be used to detect singularities of the 
Lorentzian metric of the Universe using measurements of the Cos¬ 
mic Microwave Background radiation. The physical model we 
study is mathematically rigorous but highly idealized. 


1. Introduction 

We study the dectection of singularities of the Lorenzian metric of 
the Universe from Cosmic Microwave Background (CMB) radiation 
measurements. The singularities are considered in the sense of the wave 
front set that describes where the metric is non-smooth in the spacetime 
and also in which direction the singularity occurs. The direction of the 
singularity is characterized by using the Fourier transform of the metric, 
see Definition 12.11 below. 

A singularity in the metric could be caused for example by a cosmic 
string 0 [33]. A cosmic string is a singularity in the stress energy 
tensor that is supported on a two-dimensional timelike surface in the 
spacetime. The existence of cosmic strings finds support in super¬ 
string theories ISj, however, there is no direct connection between 
string theory and the theory of cosmic strings. We refer to mum 
regarding the existence (or inexistence) of cosmic strings in view of 
CMB measurements collected by the Planck Surveyor mission in 2013. 

The singularities of which potential detectability is interesting to 
study include cosmic stings, monopoles, cosmic walls and black holes. 
There is a vast physical literature concerning the effects of particular 
types of singularities or topological defects on the CMB measurements, 
see e.g. mum and references therein. The contribution of the present 
paper is to adapt techniques from the mathematical study of inverse 
problems to CMB measurements. These techniques allow us to detect 
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singularities without apriori knowledge of their geometry. Hence it 
might be also possible to detect singularities that are not predicted 
by the current physical knowledge. Furthermore, the techniques allow 
us to study the opposite question, that is, what type of singularities 
are invisible to our measurements and therefore can not be detected 

mm- 

Several types of measurements have been proposed in astrophysical 
literature for detection of topological defects of Universe. These include 
optical measurements, such as gravitational lensing effects caused by 
cosmic stings ®mm, observations of density of mass in Universe, 
measurements of gravitational waves |ljj, and the temperature changes 
in the Cosmic Microwave Background mmmm- We show that the 
CMB measurements have a tomographic nature, and concentrate on 
them in this paper. 

The detection of singularities has been extensively studied in mi¬ 
crolocal theory of tomography, see e.g. [12, 132, 22] and the review 0. 
In many problems related to tomography, for instance in medical imag¬ 
ing, it has been shown that measurements can be used to detect sin¬ 
gularities. Moreover, the visible singularities have been characterised 
in many cases. In this microlocal context, a singularity is considered 
to be invisible if it causes only a smooth perturbation in the measured 
data. 

In the present paper we characterize the singularities of the Lorenzian 
metric that are visible from a linearization of the CMB measurements, 
and that move slower than the speed of light. We show also that all 
singularities moving faster than the speed of light are invisible. We 
do not analyze recovery of singularities moving at the speed of light. 
Our approach is based on a highly idealized deterministic model of 
CMB measurements, but such a model can be viewed as a first step in 
developing tomographic methods for more realistic, possible stochastic, 
models of the CMB measurements, see e.g. mi- 

We obtain the characterization of visible singularities via microlo¬ 
cal analysis of the geodesic ray transform on a Friedmann-Lemaitre- 
Robertson-Walker type spacetime. The transform is restricted on light 
rays and we call it the light ray transform. Such a light ray transform 
has been studied in 1 + 2 dimensions in [16] , and the 1 + 3 dimensional 
light ray transform, as considered in the present paper, belongs to the 
class of Fourier integral operators studied in O' 

Contrary to [IT], we avoid using the I P:l calculus by first microlo- 
calizing on spacclike covectors. Then we invert microlocally the light 
ray transform up to potential fields, see 0 below for the definition, 
and conformal multiples of the metric of the spacetime. This is sharp 
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since those two subspaces belong to the kernel of the linearization and 
they correspond to the gauge invariance of the non-linear problem un¬ 
der diffemorphisms and conformal changes. The spacelikc covectors 
correspond physically to singularities moving slower than the speed of 
light. 

This paper is organized as follows. Section 1 is the introduction 
and Section 2 introduces some notations. In Section 3 we formulate 
the inverse problem for the CMB measurements and its linearization, 
and in Section 4 we state the main results. Section 5 deals with the 
parametrization of the CMB measurements. In Section 6 we describe 
the conformal invariance inherent in the problem. Section 7 contains 
the reduction of the linearized problem to inversion of the light ray 
transform. In Section 8 we study the light ray trasform in a translation 
invariant case and express its normal operator as a Fourier multiplier. 
This motivates our subsequent study in the general case which is not 
translation invariant due to the fact that measurements are available 
only in a small set in the spacetime. In Section 9 we study the null space 
of the light ray transform that corresponds to the gauge invariance of 
the problem. In Sections 10 and 11 we characterize the visible spacelike 
singularities. Moreover, we compute the symbol of the normal operator 
of the light ray transform on the cone of spacelike covectors, on which 
it is a pseudodifferential operator. 


2. Notations 

Let M C R n be open. We denote by V (M) the distributions on 
M and by S'(M) the distributions with compact support. Moreover, 
we denote by C m (M), m = 0,1, 2,..., oo, the space of m- times con¬ 
tinuously differentiable functions, and let C™(M) C C m (M ) be the 
subspace of functions with compact support. For m < oo, we define 
also the Banach space 

C™(M) = {ue C m (M ); V sup |<9>(x)| < oo}. 

\a\<m XGM 

We denote by Sym 2 the symmetric 2-tensors on M, and by A 1 the 
1-forms on M. We use the notations F(M\E), T = V, S', C m , C™; 
E = Sym 2 , A 1 , C, for the subspaces T of distributions taking values on 
the vector bundle E. 

We denote by S m (T*M ) = S™ 0 (T*M), m G M, the symbols as de¬ 
fined in [T8J. If y G S m (T*M ) then we denote the corresponding pseu¬ 
dodifferential operator also by x- We denote the wave front set of a 
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distribution u G T>'(M]E) by WF(tt). For reader’s convenience we 
recall the definition here. 

Definition 2.1. Let u G V'(M] E), E = Sym 2 , A 1 , C. A point (x, £) G 
T*M \ 0 is m the wave front set of u if there does not exists an open 
cone r C TfM containing £ and a function 0 G Cq°(M) satisfying 
0(x) 7 ^ 0 such that for all N = 1,2,... there is C > 0 such that the 
Fourier transform of (j)u decays as follows: 

|S(OI <G (1 + |^|)- W , eer. 

Let us recall a typical example, see e.g. j20]. Th. 8.1.5]. Let u G 
C“(M), k < n, and let us consider coordinates (x',x") G M fcx ( n-fc ) 
on M. Then the Euclidean surface measure S(x") of {(x'jO^x' G 
satisfies 

WF(u<S(x")) = {(x',0;0,e") G T*M \ 0; (x',0) G supp(u)}, 

where £" is the covector corresponding to x". As all smooth surfaces 
in M are locally of this form, it can be seen that the wave front set 
of a surface measure is the conormal bundle of the corresponding sur¬ 
face. This applies in particular to cosmic strings that are singularities 
supported on two dimensional surfaces, see e.g. [53, Chapter 6 ]. 

3. Mathematical formulation of the CMB measurements 

3.1. The data and the inverse problem. Let us begin by formu¬ 
lating a mathematical problem on a time-oriented Lorentzian manifold 
(M, g) of dimension 1 + n. Below, we will consider a linearized version 
of this problem, and solve the linearized problem only in a microlo¬ 
cal sense. We assume that M is a smooth manifold and that g is a 
C 2 -smooth metric tensor with signature (—, +,..., +). 

We recall that a point (x, V) on the tangent bundle TM is called 
spacelike if the inner product (V, V) g is strictly positive, timelike if 
(V, V) g < 0, and lightlike if (V, V) g = 0. A submanifold E C M is 
called spacelike if all the tangent vectors in TE are spacelike, and a 
geodesic 7 : [0,1?] —> M is called a null geodesic if its tangent vector 
7 (r) is lightlike for one, and hence for all r G [0, P\. Let us also recall 
that (x,Z) G TM is called an observer if Z is future pointing and 
(Z,Z)g — — 1 . 

Let E C M be a smooth spacelike submanifold of codimension one, 
and let v be the future pointing normal vector field on E that is of 
unit length in the sense that (V, u) g = —1. Let U be another smooth 
submanifold of M and let Z be a smooth section of TM defined on L/, 
and suppose that (x, Z) is an observer for all xGW. 
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Let E o >0 and let us consider a null geodesic f3 satisfying 

(1) P(0) G E, 0(O) : u) g = -E o , 
and suppose that j3(r) GW for some rGi. We define 

(2) x = 0(t), V = E~ 1 /3(t)-Z, 


where E = —(/3(t),Z) 9 . Note that V satisfies 

(3) (V,Z), = 0, (V, V) a = 1, 

and that (x,V) determines uniquely a null geodesic (3 satisfying ([Tj). 
Indeed, (3(r) = x and /3(r) = E(V + Z) determine a null geodesic up 
to an affine reparametrization, that is, up to a choice of r and E, and 
([Tj) fixes the parametrization. 

We denote by S>g iX ,z the set of vectors V G T X M satisfying (J3]) . The 
set §>g tX ,z is called the celestial sphere of the observer (x,Z), see e.g. 
[21] • Furthermore, we denote by S 9i zU the set of points (x,V) G TM 
such that x G 77, V G §> g , x ,z, and that there is a null geodesic /3 
satisfying ([TJ) and ([2]) for some r Gl. 

Physically, the null geodesics (3 satisfying (JTj) correspond to photons 
that are emitted with a fixed energy E 0 uniformly in all future pointing 
light like directions on E. Moreover, E and V in ([2]) are the energy and 
the Newtonian velocity of f3 as measured by the observer (x,Z). The 
proportional difference between the emitted and observed energies, 


(4) 7 Z g>z (x,V) = 


E 0(0), v) 


I; (x, V) ^ Sg,zM, 


E Wt) ,Z), 
is called the redshift of (3 as measured by (x,Z). ffere f3 is the null 
geodesic satisfying ([Tj) and ([2]) for some rGl. A general formulation 
of the inverse problem that we consider is the following: 


Problem 3.1 (Inverse problem for redshift measurements). Given the 
function 7 Z 9t z : S g zM -G M determine (W, g ) where W C M is the 
union of the null geodesics connecting points ofU to points o/E. 


The red shift measurements (JTJ) do not change under the group of 
transformations generated by diffeomorphisms and conformal changes. 
Indeed, if if is a diffeomorpism that is identity on E and U, then the 
redshift data corresponding to g and its pullback if*g are the same. 
Next, if c is a smooth strictly positive function on M satisfying c = 1 
on E and 77, then the redshift data corresponding to eg and g are the 
same. This is based on the fact that the null geodesics of eg are the 
same as those of g as point sets but they are parameterized differently, 
see Lemma 6.1 below. Hence the determination of ( M,g ) in Problem 
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3.1 should be understood modulo the gauge invariance given by this 


group of transformations. 


3.2. The linearized CMB inverse problem. Let us now turn to 
a linearization of Problem 13.11 at a Friedmann-Lemaitre-Robertson- 
Walker type model, that is, at a Lorenzian manifold of the following 
warped product form 

(5) g(x ) = — dt 2 + a 2 (t)dy 2 , x = (t, y) G M = (0, oo) x M 3 , 

where a > 0 is smooth on (0, oo) and may have singularity as t —> 0. 
A singularity at the boundary t = 0 corresponds physically to the Big 
Bang. In particular, the choice a(t) = f 2//;! gives the Einstein-de Sitter 
cosmological model [MJ p. 31]. Occasionally we write a(x) = a(t ) for 
x = ( t,y). 

Let to > 0 be small and consider the surface 

(6) £ = {(t,y) G M : t = t 0 }. 

Then the photons satisfying (JTj) give a model for the CMB radiation. 
Physically, the time to corresponds to the time when the Universe was 
cool enough so that stable atoms could form. As these atoms could 
no longer absorb the thermal radiation, the photons that were around 
stayed forming the CMB radiation. 

The model is highly idealized as the initial frequency spectrum of 
microwave background radiation would more likely follow e.g. Planck 
photon distribution function, see e.g. m Ex. 5.5.4], A physically 
more realistic model would allow the photons emitted at the surface 
£ to have an energy distribution with median Eq. This would lead 
to similar considerations as below, by assuming that we measure the 
energy distribution of photons and find the median of this distribution. 
However, for simplicity we consider the case when all emitted photons 
have constant energy. 

Let g e G C 2 (M; Sym 2 ), e G [0,1], be a one parameter family of 
Lorentzian metrics, and suppose that g e = g in M \ Mo where 

Mq = (f 0 , oo) x M 3 . 

Here to is as in (|6]) . We will now describe the CMB measurements on 
(M,g e ). Let U > to, let Mi C M 3 be open and bounded, and define 
U = {ti} x U x . To avoid technicalities in the exposition, we assume 
that g e = g on U and consider the observers (x,d t ), x G U. We will 
study more general observers in Section [5} 
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Figure 1. The observer ( y,dt ) measures the energy of 
the photon /3. The photon has the Newtonian veloc¬ 
ity v on the celestial sphere of (y,d t ). We show that 
the measurement carries information on the singular¬ 
ity (x,f) G WF(<9 e g e | e=0 ) tha t satisfies /3(ti) = x and 
£(/3(ti)) = 0, see Theorem 


4.4 


Let S 2 be the Euclidean unit sphere in 1 
pliism 

(7) 9 : S 2 ->• S geiXidt ~ d t , 9(v) = -(1, a(fi)~ 1 u) G 


We define the diffeomor- 


D 1 + 3 


If g e is sufficiently close to the background metric g then Lemma 3.2 
below says that S g ^U can be parametrized by U\ x S 2 , and the red- 
shifts Tlg^o t of CMB photons as measured by (x,d t ), x G U, can be 
parametrized as follows 


(8) R e (y,v) = (%(T e (x,9)-,x,9),d t ) ge - 1, (y, v) G U x x S 2 , 


where x = (ti,y), 0 = 9(v), 7 e (-;x,0) is the geodesic of ( M,g e ) with 
the initial data (x,9), and 


(9) t £ (x,9) = min{r > 0; 7 £ (t;x,9) G S}. 

Lemma 3.2. Let g be a Lorentzian metric tensor of the form 0- 
Let t\ > to, let U\ C M 3 be bounded, and define U — {fi} x li\. Let 
g e G C 2 (M; Sym 2 ), e G [0,1] 7 be a one parameter family of Lorentzian 
metrics satisfying g 0 = g. Suppose that the map e (->• g e is differentiable 
from [0,1] to C 2 (M 0 ; Sym 2 ) and that g e = g in U U E. Let (x,v) G 
U x S 2 , and define 9 = 9{v) and r e = r e (x,9) for small e. Let c > 0. 
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Then the geodesic 

Pe{p) = 7 e(-cp + r e ;x,9), 

satisfies (3 e ( 0) G S and /3 e (r e /c) = x. In particular, if we choose the 
parametrization c = E 0 /(%(T e ',x,9),dt) ge , then (/3 e (0), d t ) ge = —E 0 and 
the redshift of f3 e as measured by (x, d t ) satisfies 

Kgc,d t ( x ,d+ d t ) = ( 7 e (r e -,x,9),d t )g e - 1. 


Pro of. We will show that r e is well-defined for small e > 0 in Lemma 
below. We have /3 e (0) G E, f3 e {r e /c) = x, (f3 e (0), d t ) 9e = —E 0 , and 


5.1 


the redshift is 

(&(0 ),d t ) 9e 


(P e {T e /c),d t ) gt 
since (9, d t ) 9e = 1. 


1 = 


-En 


-c(9, d t )g e 


1 = ( 7 e (T e ;x,9),d t )g € - 1, 


□ 


We are ready to formulate the linearized version of Problem 3.1 that 
we will consider. 


Problem 3.3 (Linearized microlocal CMB inverse problem). Given 
d e R e (x,v) | e= o, ( x,v ) G UxS 2 , determine the wave front set WF(<9 e g e | e= o) 
on T*W where W C M is as in Problem \3.1[ 

We characterize the spacelike vectors in WF(<9 e g e | e=0 ) that can be 
recovered, and show that all timelike vectors are smoothed out. Phys¬ 
ically, the timelike vectors in WF(<9 e g e | e= o) correspond to singularities 
moving faster than light, whence we do not expect such singularities 
to be present to begin with. 

Lightlike singularities could correspond to gravitational waves. Their 
recovery from the anisotropies in the Cosmic Microwave Background 
radiation is an interesting question. However, we leave this as a topic 
for future work. 


4. Statement of the results 

Let / be a symmetric 2-tensor on a Lorentzian manifold ( M,g ), and 
define the light ray transform of / by 

(10) X g f{x,9)= [ ( fi m °l)i l i m dr , 

Jr 

where x G M, 9 G T X M is lightlike, and y(r) = 7(77 x, 9) is the geodesic 
of ( M,g ) with the initial data (x,9). When g is of the form (|5]) , we 
define for t\ > to, 

L tl f(y, v) = Xgffx , 9), ye M 3 , v G S 2 , 
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where x = (h,y) and 9 = 9(v) is defined by ({7]). Note that for an 
arbitrary Lorentzian manifold (. M,g ), the map X g may fail to be well 
defined even when / G Sym 2 ). This is the case, for example, 

if (M, g ) has closed null geodesics. However, we show the following 
lemma. 


Lemma 4.1. The light ray transform L tl is a Fourier integral operator 
of order —3/4 whose associated canonical transformation is the twisted 
conormal bundle of the following point-line relation 

{( p,y,v) G M x M 3 x S' 2 ; p = 7 (r; (ti,y),9(v)) for some t G M}. 

In particular, L tl is continuous from S'(M; Sym 2 ) to 1T(M 3 x S 2 ). 


The proof of the lemma will be presented in Section 10 The following 


theorem reduces the linearized CMB inverse problem to inversion of a 
limited angle restriction of the light ray transform. 


Theorem 4.2. Let t\ > t 0 , let li\ C M 3 , and let g e G C 2 (M; Sym 2 ), 
e G [0,1], satisfy the assumptions of Lemma\3.^\ Define 


f{t,x') = 
Then on U\ x S 2 , 


d e g e (t,x% =0 , te(t 0 ,ti), 
0 , otherwise , 


d e R f L =n = 


Lt 1 a £ a d t a ~ f i 


2 a(t 0 ) 


where C a d t is the Lie derivative along the scaled world velocity ad t , 
and the powers a 2 and a~ 2 of the warping factor are interpreted as 
multiplication operators. 


The light ray transform has a non-trivial null space. 

Lemma 4.3. Let g be a Lorentzian metric tensor of the form ([3Jj and 
let ti > 0. Define Mi = (0,fi) x M 3 and 

J\f = {eg + d s u ; c G £\Mi ), u G £'{M \; A 1 )}, 


where d s is the symmetric differential defined in coordinates as follows 


( 11 ) 




(\7 jUj)j -(- (X jUj)i 
2 


Here V, = Vq t is the covariant derivative on (M,g). Then L tl f = 0 
for all f G Af. 
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Proof. By density of Cq° in £', we may assume that c G and 

that u j G Cq°(Mi; A 1 ), and consider f = eg + d s LO. Then L tl f = 0 on 
M 3 x S 2 , since ( 7 , 7)3 = 0 holds for null geodesics 7 and 

d T uj(j(r)) = (d s o;)y7 l (r)7 J (r) 

holds for all geodesics 7 . □ 


In fact, this lemma holds for every Lorentzian metric as far as the 
null-geodesics are non-trapping. Note that the group of transforma¬ 
tions g i-A- g + h, h G AT, is the linearization of the gauge invariance of 


the nonlinear problem, see the discussion after Problem 3.1 


The lemma is an analog of the corresponding results for compact 
Riemannian manifolds, where it is known that the geodesic ray trans¬ 
form vanishes on potential fields, i.e., on tensors of order m > 1 which 
are symmetric differentials of tensor fields of order m — 1 vanishing on 
the boundary, see [29J. What is new here is the scalar multiple eg of 
the metric. 

Our main result is the following. 


Theorem 4.4. Let g be a Lorentzian metric tensor of the form ©• 
Let t,\ > 0, let U\ C M 3 be open, and define Mi = (0, ti) x M 3 . Let 
(x,f) G T*Mi be spacelike, and suppose that there is a null geodesic 7 
of (. M , g) and Ti, 72 G R such that 

(12) 7(o) = x, £(7(ti)) = 0, 7 ( 7 - 2 ) G {ti} x Ui. 


Then there is x £ S°(T*(M. 3 x S 2 )) supported in T*(IA\ x S 2 ) such that 
for all f G £'(M] Sym 2 ) the following are equivalent 
( 1 ) (*,£)GWF (L* tiX L t J), 

(ii) (x,£) G WF(/ + h) for all h G M■ 


In particular, the theorem says that the operator Lf xLt, does not 
move spacelike singularities. We will give an explicit choice of x hr 
Section 


11 


and show that L* ti \L tl is a pseudodifferential operator of 
order — LvVe will also compute its principal symbol after a conformal 
scaling, see Proposition 11.4 below. 

Note that the visibility condition (12) is analogous with the visibility 

Tlr. 3.1]. 


condition for limited angle X-ray tomography, see e.g. |22J Th. 3.1]. It 
is also sharp in the sense that if for a spacelike (x, £) G T*M there does 
not exist a null geodesic 7 satisfying ( 12 ) then for all x £ C^ilti) the 
wave front set WF(yL il /) is empty if WF(/) is contained in a small 
conical neighborhood of (x, £) in T*M. 

We will also show that all timclike singularities are lost in the sense 
that if WF(/) contains only timclike covectors then WF(L tl /) is empty. 
In other words, L tl is smoothing on the cone of timclike covectors. 
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5. Parametrization of the CMB measurements 


In this section we show that the function r e , see |f]), is well-defined. 
We do this in a slightly more general context than that in Lemma 
2 , that is, we relax the assumption that g e = g on U, and consider 
observers (x, Z e ) on U where Z t is not necessarily d t . 

Let g e G C 2 (M; Sym 2 ), e G [0,1], be a one parameter family of 
Lorentzian metrics satisfying g 0 = g, where g is a Lorentzian metric 
tensor of the form ([5]). Let W e C M be a spacelike submanifold of 
codimension one, and let Z e be a section of TM defined on W e such 
that (x, Z e ) is an observer with respect to g e . Let U\ C M 3 be open 
and suppose that <f> e : U\ —* W e is a diffeomorphism. In this section we 
consider measurements by the observers 


(13) 


M’..(?/). 2,). yeUi. 


We assume that Wo = {ti} x U\ where t\ > to, $o (y) = and 

that Zq — d t . 

For small e > 0, the CMB measurements by observers (x, Z e ) can be 
parametrized by the vectors in the celestical spheres § g e , x ,z e , x G W e . 
Indeed, we have the following lemma: 


Lemma 5.1. Let M C M n be open, let g be a smooth Lorentzian metric 
tensor on M, let Z be a smooth vector field on M, and suppose that 
(x, Z) is an observer for all x G M. Let Kq C M be compact and define 

K = {(x, V-Z)e TM- x G K 0 , V G 

Let Ec M be a smooth submanifold of codimension one, and suppose 
that the geodesics 7 (-;x, 6 l ), (x, 9) G K, intersect £ non-tangentially. 
Then there are neighborhoods Q C C 2 (M; Sym 2 ) of g and K, C TM of 
K such that the geodesics 7 h(-',x,9), (x,0) G /C, with respect to h G Q 
intersects*. Moreover, the function 

T h (x, 9) = min{r > 0; 7 h (T;x,d) G £} 

is C 1 in /C x Q, and the functions 7 ^(r; x, 9) and 7 ^(r; x, 9) are C 1 in 
a neighborhood of the set 

{( r,x,9,h ); r G [0, t^(x, 9)\, (x, 6 i )g/C, heQ}. 

Proof. Let Q 0 C C 2 (M] Sym 2 ) be a neighborhood of g such that all 
h G Go are Lorentzian. We define the open set U = TM x Go 011 th e 
Banach space E = TM x C 2 (M; Sym 2 ), and consider the vector field 
F : U -7 E, F(x, 9, h) = (9, /(x, 9, h), 0), where / : E -7 M n has 
components 

ffix, 9, h) = -r^(x, h)9 k 9 l , (x, 9) G TM, h G C 2 (M; Sym 2 ), 
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and r^(x, h ) are the Christoffel symbols of the metric tensor h at x. 
We recall that 




1 j^jm f dfornk dflmg dhkg\ 

2 y dx k dx m ) ’ 


where W m is the inverse of hjk, and see that Y 3 H : M x Qq —> M is 
continuously differentiable. Hence F is continuously differentiable, and 
the flow of F is continuously differentiable on its domain of definition 
Dorn(F) C lx [/, see e.g. [21 ( Th. 6.5.2], Note that the projection 
of the flow on TM gives the geodesic flow with respect to the metric 
tensor h. 

Let k = (x,9) G K. We have assumed that (t 9 (a), k, g) G Dorn(F). 
As 7 (-;k) intersects S non-tangentially, the implicit function theorem 
gives neighborhoods /C K C TM of n and Q K C C%(M] Sym 2 ) of g such 
that Th is C 1 on fC K x Q K . The set K is compact since K 0 is compact and 
S g , x ,z is diffeomorphic to the sphere S n ~ 2 . We choose a finite set J C K 
such that /C K , k G J, is an open cover of K , and define /C = U Ke j/C K 
and ^ = n K£j Q K . □ 


Let us now continue our study of the observers (13). We choose a 
smooth family of diffeomorphisms 'L x ,c : S 2 — > S 9 e , 3 ;,z, £ G W e , such 
that 

= (°> -a(F) _1 v) e K 1+3 , x G W 0 . 

Then the redshift measurements can be parametrized as 


R*(y, v) = ( 7 e(r e (x e , 0 e ); x e , 0 £ ), Z £ ) ge , (y, v) G U x x S 2 , 

where x t = & e (y), 9 e {x,v) = ^ x<iie {v) — Z e , 7 £ (-;x, 0) is the geodesic of 
(M,g e ) with the initial data (ay 0), and r e is dehned by (|9]) . Note that 
9 0 = 9 where 9 is dehned by ([7]). 

We omit writing e as a subscript when it is zero. The linearized 
measurements satisfy 


d e Re(y, , y)| e =o = d e R e (y, n)| e=0 + d e ('y(r(x e , 9 e )] x e , 9 e ), Z e ) g | e=0 , 

where the second term is a smooth function on U\ x S 2 if the derivatives 
at e = 0 of Z e and the local parametrizations and e are smooth on 
IA\ x S 2 . Thus the second term plays no role in the microlocal analysis 
of d e R e (x,v) |e=o an d we are reduced to the case covered in Theorem 

m 


6. Conformal invariance 

In this section we show that Theorem I4.4l is invariant under conformal 
scaling. Let us recall that a metric tensor g of the form (J 5 J) is conformal 
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to the Minkowski metric tensor in suitable coordinates. Indeed, we 
define the strictly increasing function 

(14) s(t) := f a(t) _1 (if. 

Jo 

Then 

g(s, x') = a 2 (s)(—ds 2 + ( dx ') 2 ), (s, x') G ( 0 , oo) x M 3 . 

For example, the Einstein-de Sitter model corresponds to aft) = t 2 / 3 , 
and then s(t) = ?>t 1 ^ and a(s) = s 2 / 9. 

It is well known that null geodesics are invariant up to a reparametriza- 
tion under a conformal change of the metric. For the reader’s conve¬ 
nience we give a proof here. 


Lemma 6.1. Let g be a smooth Lorentzian metric tensor, let a be 
a smooth strictly positive function, and define g = a~ 2 g. Then the 
null geodesics 7 of g correspond to the null geodesics p of g under the 
reparametrization 

( 15 ) 7(r) = /i(cr(cr)), r(o-) = / a(p(p)) 2 dp , 

Jo 

where c = a(7(0)) 2 , r(cr) is a strictly increasing function and cr(r) is its 
inverse function. The reparametrization is chosen so that 7(0) = p( 0 ) 
and 7(0) = /i(0). 


Proof. We define a by e a = a. Koszul formula implies 


(16) V = V A W + (Xa)Y + ( Ya)X - (X, Y)~ g , grad^a, 

where grad^ is the gradient with respect to the metric g. We apply this 
formula to X — Y — p, where p is a null geodesic with respect to g. 
Then the first and the last term on the right hand side vanish and we 
have 

D a p = 2 (a o p)'p, 

where D a p is the covariant derivative of p with respect to g along p , 
and the prime denotes the derivative of a real valued function on M. 
We have r' = e 2 " 0/i , o' = e ~ 2ao f 1 ° a an( i 


a 


= — 2(a o p o a)'a' = — 2 ((a o p)' o a)(a')‘ 


Let us consider the curve 7 defined by the formula (15) with c = 1. 
Then the covariant derivative D r 7 with respect to g along 7 satisfies 


D r 7 = D T (a\p o a)) 

= d T (a'(p k o a))d k + {cr') 2 {p l o a)(p J o cr)(r^ o p o a)d k 
= a"(p o a) + ( a') 2 (D a p ) o a = 0. 
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Here T*- are the Christoffel symbols of g. We see that 7 is a geodesic 
with respect to the metric g. The same is true for any 0 since affine 
reparametrizations of geodesics are geodesics. Moreover, 

(17) 7 (t) = ca“ 2 ( 7 (r))/i(p), p = a(cr). 

□ 

Corollary 6.2. Let g be a smooth Lorentzian metric tensor, let a be 
a smooth strictly positive function, and define g = a~ 2 g. Then the 
corresponding light ray transforms satisfy X„f = a 2 X„a^ 2 f. 


Proof. We set p = a(cr) and use (17) and dr = c l a 2 dp. □ 

Lemma 6.3. Let g be a smooth Lorentzian metric tensor, let a be 
a smooth strictly positive function, and define g = a~ 2 g. Let c be a 
smooth function and uj be a sm.ooth 1-tensor. Then 

a~ 2 (cg + d s u> ) = eg + d s uj , 

where c — c + (uj, d\oga)g and uj = a~ 2 u). Here d s is the symmetric 
differential with respect to g. 


Proof. We define a by e a = a, and write (|16|) in coordinates. This gives 
the following fc 
g, respectively, 


the following formula for the Christoffel symbols I 'b and I 'b of g and 


r£ = r f f j + sfd j a + s , ;d l a-b k g lj , 

where b k = (gradual) fc . The symmetric derivative transforms as 

(d s u))ij = ^ (diUJj + djuf) - T \jUJk = ( d s uj)ij - - ujjdia + uj k b k g i:j . 

Moreover, 


e 2 a (d s (e- 2 a o;))y = (d s cofi 


Ujdia 


Uidja. 


□ 


Let g be of the form (|5]) and define g = a~ 2 g. Let x G M and let 
£ G T*M be spacelike for g. Then £ is spacelike for g. Suppose that 
7 is a null geodesic of (M, g) and that there are 0 < 17 < r 2 such that 
(12) holds. Let p. be the reparametrization (15) of 7. Then p is a null 
geodesic of ( M,g ) and there are pi,p 2 GM. such that 

P(pi) = x, £(/t(pi)) = 0 , 7 (p 2 ) e {ti} x Ui. 

Suppose now that Theorem 4.4| holds for 'g, that is, there is a sym¬ 
bol of order zero x supported in T*(14\ x S 2 ) such that for all / e 
£'(M] Sym 2 ) the following holds: 

(r,£)6WF(L* X Z tl /) iff (x,£)GWF (f + cg + d 8 u) 
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for all c G V(M) and oj G V(M; A 1 ). Here L tl (y, v) = Xg(x, 9 ), where 
x = (ti,y) and 9 = 9(v) is defined by ([7). Note that L tl = a 2 (ti)L tl a 2 , 
and 

L * tl X L ti = a^iti)^ 2 L* ti xL tl a~ 2 . 

Let / G S'(M-,Sym 2 ), c G V\M) and c o G V{M\ A 1 ), and let us 
define / = a -2 /, u) = a~ 2 uj, and c = c + (uj, d\oga)g. Then 

cl T cy T d ui) — f T cy T d uj y 

and 


(x,0 G WF(L* tiX L tl f), 
(x, 0 G WF(/ + c(7 + d s (n) 


iff 

iff (s,()GWF(/ + ^ + d s w). 


Thus Theorem 4.4 holds for g if it holds for g. 


7. Derivation of the tomography problem 


In this section we prove Theorem 4.2 


The proof is similar to that 
in [23]. All the differentiations below are justified by Lemma 5.1 We 
omit writing e = 0 as a subscript. 

Let t\ > t 0 , let U.\ C M 3 , and define U = {ti} xU\. Let (x, v ) ElAxS 2 
and write 6 = 9(x,v). Note that d t = a~ l d s where s is dehned by (14). 
We have 


d e {i e {r e (x,vy,x,9),d t )g e — a 1 (t 0 )d e (%(T(x,v)]x,6),d s ) ge , 
since 7 = 0 and therefore 

d e ( X (r e (x, n); x, 6), d s ) g = d T ( 7 (r; x, 9), d s ) g \ T=Te ^ v )d e T e (x, v) 

= (j(T e (x,v)-,x,9),d s ) g d e T e (x,v) = 0. 

Moreover, 

a e (7 e (0; x , 9), d s ) ge = d e (9, d s ) g = 0, 
since g e = g in U. Hence 


(18) a(t 0 )d e R(x, v ) = d e ((7 c (t(x, v);x, 9), d s ) 9t - ( 7 e ( 0 ; x, 9), d s ) 


9e) 


We use the reparametrization (15) and recall the formula (17). Then 
(19) 


=c(A«(P \x,0),d, 


s)g t i 


where c = a(x) 2 , p = cr(cr), g e = a~ 2 g e and y e (p;x,9) is the geodesic 
of ( M,g e ) with initial data (x,9). 

The computations below will use the fact that g = —ds 2 + (dx') 2 
is a constant tensor in the (s, x') coordinates. When computing in 
coordinates we write (s,x') = (x°,x 1 ,x 2 , x 3 ), and use also the notation 
de{p) = Pe(p]x,9). 
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The geodesic equations for p, e can be written in the form 


( 20 ) 


1 


9p(g £ JkP £ ) — 2 ipxi9e,lm) ft> £ ft> £ 


indeed, the equation (20) can be transformed to the equation, 

fyi + rijJ&T = o, 

written in terms of the Christoffel symbols T J elm of the metric g £ , by us¬ 
ing the chain rule on the left hand side and symmetrizing the obtained 
sum, see e.g. [32] p. 51]. 

We denote h = d £ g £ | e=0 and have 

d e {(d xJ g eM ) | e=0 = (. d x jd e g eM \ e=0 ) p}fi m + ( d x ig im ) d £ (/4/i™) | e=0 

= (■ d x ih lm )fi l fi m , 

where we have used the fact that g is a constant tensor. Hence 

9pd £ (fl £j ^ ) g, c 0 9 e 0p (g £ , Ofc/b) |e=0 2 (^x°^lm) /I /i 

We integrate this with respect to p and obtain 

(21) ^ ( d s hi m )p l fi m dp = d e (fi £ (r),d s )g e \ e=0 - d £ (fi £ (0),d s )g e \ e=0 , 

where r > 0. Now (21), (18) and (Jl9|) imply 


c 1 a(t 0 )d e R(x,v) | £=0 = - / ( d s h lm )fi l p m dp , 

^ Jo 

where cr 0 = <j(ct(x, v)). 

Let t\ > 0, let / be the cutoff of d £ g £ | e= o as in Theorem 4.2, and 
suppose that x € Lb. Then h = a~ 2 f on /i([0,cr 0 ]) and / = 0 on 
g(R \ [0, cr 0 ]). Note also that (C a a t h)i m = (C da h)i m = d s hi m . Hence 

d e R(x,v )| e=0 = n _ C , XgC adt a~ 2 f = ] X g a 2 C adt a~ 2 f, 


2 a(t 0 ) 

and we have proven Theorem |4.2 


2a (to)' 


8. Backprojection of the light ray transform 

Let us now turn to the proof of Theorem |4.4 By the conformal in¬ 
variance, we may assume that the background metric is the Minkowski 
metric, that is, a = 1 identically. If we have data on all the light 
rays, i.e., if IA\ = M 3 , then the light ray transform L tl is invariant with 
respect to translations. Although the case U\ = M 3 is physically un¬ 
realistic, we discuss it briefly in this section since it allows for very 
explicit computations. 
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After a translation in the t coordinate, we may suppose that t\ = 0. 
We write L = Lq, and have 

(22) Lf(y,v) — [ fim(s,y + sv)6 l 6 m ds, y G M 3 , v G S 2 , 

Jr 


where 9 = (l,u) € R 1+3 . Note that we have made change of variable 
s = —t in comparison with (10), and also use 9 = —9{v) in comparison 
with (Jt|) . 

We use the Euclidean surface measures in all the integrations below. 
Let f,he C^(R 4 :Sym 2 ). Then 


(■ Lf , Lh) L 2 ( R 3 xS 2 ) 



is 2 


fi m (r, y + rv)hjk(s, y + sv)9 4 9 k 9 l 9 m dsdydrdv 



is 2 


fim(r, x ' + ( r — x°)v)hjk(x)9 2 6 k 9 l 9' n dxdrdv, 


where we have used the fact that 


(23) (s, y) i—> (s, y + sv) — x — (x 0 , x') 

is a linear isometry on M 1+3 . After making the change of variables 
p = x° — r, we see that and the normal operator L*L is the convolution 
K^ klm * firm with the kernel 


(I< 4klm ,cP) v , xvm = I [ 9 j 9 k 9 l 9 m (f>(p9)dpdv. 

Js 2 Jr 

The following lemma is in the spirit of the representation of the ray 
transform in [29l [31 ] if we allow for a singular weight on the light cone. 

Lemma 8.1. The Fourier transform of K ]klni is a locally integrable 
function and satisfies 


ffjklm \ _ 


2 tt (|^' / | 2 — |<Co| “) 1 ^ 2 f s i 9 j 9 k 9 l 9 m dv, f is spacelike, 
0 , otherwise. 

Heref=(f 0 ,f')eR 1+3 ,9 = (l,v), | • | is the Euclidean norm, and 

(24) Sl = {veS 2 -, f o + f'v = 0} 

is a circle of radius |£ , | _1 '\/l£ , l 2 — |£o| 2 - 
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Proof. To simplify the notation, we fix the indices j, k, l,m and denote 
o(i?) = 6 3 6 k 6 l 9 m . Let us also define k(v,£) = £ 0 + £'v. Then 

(jW w "\0) 5 , x5(ffi4) = f [ [ ae-WWWdpdv 

J S 2 J K JR 4 

= 2vr / {k*S,4>)s'xS( r 4 ) adv. 

Js 2 

Here the pullback is by k v (£) = k(v,£) with hxed v. Note that the 
gradient Vic v — 6 does not vanish, and therefore the pullback is well- 
defined. 

Let us assume for a moment that supp (0) does not intersect the set 
of lightlike vectors 

{£ = (6),OeRxR 3 ; |£o| = |£U 

Then 

/ (ft*<5, 0)s'xS(R 4 ) a d y = ((ft* 8, fl)p'xD(S 2 )) 0)s'xS(R 4 ), 

Js 2 

Here the second pullback is by ft %(v) = /c(u,£) with hxed £. Let us 
now show that this is wcll-dehned. We have Vr 31C£ = £'. Thus on 

s l = ft^(°) 

|V 5 2 ft^(u)| 2 = |£' - (f'v)v\ 2 = W\ 2 - \?v\ 2 = in 2 - l^ol 2 . 

We see that the pullback by is wcll-dehned if and only if £ is not 

lightlike. 

If £ is timclike, then |£'| < |£ 0 | = |£ 7 u| < |£'| on which is a 
contradiction. Thus does not intersect the timelike vectors and 
Kjkim _ q on timclike vectors. 

Suppose now that £ is spacclike. Then 

(25) (ftjt<5, <x)t>'xv(s 2 ) = (|£'| 2 - |£o| 2 )- 1/2 / a{v)dv. 

The set is the intersection of S 2 with the affine plane 

Pt = {ve M 3 ; = —£ 0 }. 

As £ is spacelike, 0 < |£o| < |£'|, and we see that w = —£o£ / /|^ / | 2 G P%- 
Thus Sg is a circle of radius r satisfying r 2 + \w\ 2 — 1, that is, 

r = l ? rVl?T - l ? ol 2 - 

To hnish the proof, it is enough to show that K ]klm coincides with a 
locally integrable function. Let us hrst assume that supp (0) does not 



INVERSE PROBLEM OF FINDING COSMIC STRINGS 


19 


intersect the set 


(26) 


{£ = «„,?') eRxR 3 ; l?'l = o}. 


Then for each £ £ supp (</>) we can choose spherical coordinates so that 
v —v(a,j3 ) = (sin a cos/3, sin ctsin/3, cos a) 
and v(0,P) = — £ / /|£'|. In these coordinates 

<P(-£'v,{,')a(v)dvd£, 


is 2 


(«*<*, <j>) 


cS'x 


adv = 



's 2 J R 3 

-2?r 



R 3 JO 


0(|£'| cos a, £ / )a(a, j3) sin adadjdd^. 


'0 


The change of coordinates £ 0 = |£ 7 | cos a gives 
r r r2* r \£'\ 


IS 2 


(«*£, 0)s'xs(i> 4 ) adv = / |£‘ 


c-i 


1 0 


H?'l 


(p{£o,£')a(£o,P)d£ 0 df3d£. 


Hence away from the set (26), the Fourier transform K jklril is the func¬ 
tion 

rl-K 


(27) 


lf|- icK) / oKo 


where 1 q is the indicator function of the set of spacelike vectors, that 
is, 1(7(£) = 1 if £ is spacclike and it is zero otherwise. Note that a is 
bounded and that the function (27) is locally integrable. 

We have seen that K jklrn vanishes on timelike vectors and that it is 
a function away from the set (26). The origin is the only vector in the 
set (26) that is not timelike. As K jklm is homogeneous of degree —3, its 
Fourier transform is homogeneous of degree —1. The Fourier transform 
Kjklm a f linc ti on since the only distributions supported in the origin 
are the linear combinations of derivatives of the delta distribution, and 
these derivatives are homogeneous of degree —4 or less. □ 

Note that the integral f gl QjQ k Q l Q m dv is homogeneous of degree zero 

with respect to £. It can be evaluated explicitly for all the 4 4 combina¬ 
tions of the indices. As an example, let us give Kj J3J (^), j = 0 , 1 , 2, 3, 
for a spacelike £ £ M 4 . After a rotation in £', we may assume that 
£ = (£o, 0, 0, —1£'|). Then at £ it holds that 


A H 


0000 


47T 2 

(0 = ^ 


2\2 


A' im (£) = 
and A' 2222 (£) = A' im (£). 


37t 2 (|£T - |£o| 2 ) 

2 |£? 


AT 4ddd (£) = 


4vr 2 £ 0 4 

l£'l 5 : 
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An analogous light ray transform can be defined for scalar functions 
as follows, 

L sc f(y,v)= / f(s,y + sv)ds, y G R 3 , v e S 2 , 

Jr 

and a variation of the above argument shows that L* SC L SC is given by the 
Fourier multiplier A' 0000 . The scalar case will be considered in depth 
in [50] . 


9. The null space of the symbol of the normal operator 


In this section we will characterize the null space of the tensor 
Kjbim (£) f or spacelike £ G M 4 . Note that if to G (^(M 4 ; A 1 ), then 
the Fourier transform of its symmetric differential is 


(d a w)ij(£) 


2 


To cover also the limited angle case, to be discussed in the next two 
sections, we will prove the following slightly more general result. We 
refer to (24) for the definition of Sg. 


Lemma 9.1. Let £ G M 4 be spacelike, let \ £ C°°(S]) and suppose that 
X > 0 and that x does not vanish identically. Then the null space of 
the linear map N : fi m i-a N^ klm fi m on Sym 2 , where 


(28) N jklm = [ x o°o k e l e m dv, e = (i, v ), j, k = o, 1,2,3, 

J s\ 

is Ker(N) = {cgi m + fiuj m + £ m c up, c G M, w G M 4 }. Here g is the 
Minkowski metric. 


Proof. We have gi m 0 l d m = 0 since 0 is lightlike, and 9 m f m = 0 since 
v G Sg. Thus the tensors cg + £,<S>uj + uj®t; are in Ker(A^). Let us 
now show that there are no other tensors in Ker(A^). Suppose that 
/ G Ker(N). Then 

0= [ Xf jk 0 j e k 9 l 0 m f lm dv= [ X \fjk9 j 9 k \ 2 dv. 

Jsj Js i 

Then fjk9^9 k = 0 for v in a nonempty open subset of S^. 

We write £ = (£o,£ / ) € M 1+3 . After a rotation in £', we may assume 
that £ ; = (0,0,1). As £ is spacelike, there is a Lorentz boost B such 
that £ = Be 3 , where e 3 = (0, 0, 0,1). The boost A is a linear map B k 

such that if 9 k = B k 9i then 9 J = 9\ j = 1, 2, and 

(29) 9° = 9° cosh a + 9 3 sinli a , 9 3 = 9° sinli a + 9 3 cosh a , 
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where a is the hyperbolic angle (or rapidity) of the boost. 
When v G and 9 = (1, v) we have 

0 = = ( Be 3 )9 = 9 3 , 


since the boost B is symmetric. Combining this with (29) and 6 0 = 1 
gives 9° = 1/cosh a. Indeed, 

sinh 2 a 1 


3 smh a ~ 0 

9 — --— and 9 = cosh a 

cosh a 


cosh a cosh a 


Let us denote A = B 1 and define a 3 x 3 symmetric tensor h by 

him — cosh -2 a fi m , l,m = 0,1,2, 

where fi m = A^fj^A^. Let us also write W = cosh a(9°,9 1 ,9 2 ) G S 
and w = (W 1 , W 2 ) G M 2 . Then 


0 = f jk P9* = fi m 9 9 m = h lm W L W m , 


since 9 3 = 0. Moreover, W° = 1 and w = cosh a(0 1 ,0 2 ) satisfies 
\w\ 2 = cosh 2 a(l — |0 3 |) = cosh 2 — sinh 2 a — 1, 
since 9 = (l,u) and v G S 2 . Thus 

h((l,w),(l,w)) = h lm W l W m = 0 


for w in a nonempty open subset of S 1 . After a rotation we may assume 
that w = (sin /3, cos /3) is in this set for small /3. Differentiating three 
times with respect to /3 we get 


(30) 0 = h((l,w), (l,iu)), 

(31) 0 = h((l, w), (0, w')), 

(32) 0 = h((0, w'), (0, w')) - h(( 1, w), (0, w)), 

(33) 0 = -3h((0X)>(0,w)). 


Now (33) implies that h \2 = 0. By (33) we have 

h(( 1, w), (0, w')) = h(( 1,0), (0, w'))+h((0 , w), (0, w')) = h(( 1,0), (0, w')), 


which implies together with (31) that /io 2 = 0. Differentiating this once 
more we see that hoi = 0. As all the cross terms vanish, (30) implies 
that hn = —h 00 . Differentiating (33) we see that 

0 = h((0 , w), (0, w)) - h((0 , w'), (0, w')), 


which implies that hn = /i 2 2 - Thus h is proportional to the 1 + 2 
dimensional Minkowski metric. 
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We have shown that there is c 6 K such that f m i — cg m i vanishes 
when m,l = 0,1, 2. Moreover, we can choose u G M 4 such that 

/ — eg — e 3 <8) u — u (g> e 3 = 0. 

We recall that C, = Be 3 . As the 1 + 3 dimensional Minkowski metric g 
is invariant under the Lorentz boost B , we have 

/ = cs + f®o; + w<8>f, 

where c 0 = Buj. □ 


Let us now construct a projection onto a complement of Ker(lV) as 
a symbol P of order zero that commutes with N. The construction 
is inspired by [29, 31]. We will raise and lower indices by using the 
Minkowski metric g. We define 

x ki = 44 ; 4 = 4 - VkV\ V = t == f ; 

vy(£4) 

where £ is assumed to be spacelike, and <5[, is the Kronecker delta. Then 
rf'rj, = g(r),r]) = 1 and 

(34) K* k r}i = rj k - VkrfVi = 0. 

Thus £ $3 w + a; <g) £ G Ker(A) for any u 0 G M 4 . 

We have Pyjij = gkj — VkJlj, whence Xg — g, where 


l-iki = gki ~ VkVi- 

Let us define /W = 5 + — rfrf , and 

H®h 


(35) 


P — A — 


Lemma 9.2. Lef N be a tensor of the form (28). Then P is a projec¬ 
tion satisfying Ker(P) = Ker(lV) and PN = N = NP. 


Proof. We have 

h v 9ij = g v 9ij - g(v, v) = 4 - 1 = 3 , 

whence g G Ker(P). Moreover, 

(36) hPr] t = if — rfrfrji — 0, 


whence £ <g) u> + oj 
Ker(lV) C Ker(P). 

Let us show that P 2 = P. 

A 2 = A. Similarly (36) implies that 


£ G Ker(P) for any a; G M 4 . We have shown that 


We have 44 = 4 by (34), whence 


h lJ ga = S% - g(v, v) = 3 , 
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and (p ® h) 2 — 3/i ® h. Moreover, 

A(/i <g) h ) = (Xp) ® h = (A g) ® h — p® h. 
Analogously with above, we have p® h\ = p ® h. Thus 



We have PffNf^ = N l k ] since 



Analogously NP = N, therefore Ker(P) C Ker(N). 


□ 


Lemma 9.3. Let M C R 4 be open and let (x,C.) € T*M be spacelike. 
Then for all f G £'(M\Sym 2 ) we have that (x,£) G WF (Pf) if and 
only if (x, £) G WF(/ + eg + d s uo) for all c G £'(M ) and u> G £'(M; A 1 ). 
Furthermore, if A is a pseudodifferential operator satisfying A(cg + 
d s u) = 0 for all c G £'(M) and u G £'(M\ A 1 ), then (x, 0 G WF (APf) 
if and only if (x, f) G WF (A/). 

Proof. Let y G S°(T*M ) satisfy y(x,£) — 1 and suppose that supp (y) 
contains only spacelike covectors. Note that (x, £) G WF (Pf) if and 
only if (x, £) G WF(Py/), since P(1 — y) is smoothing near (x,£). 

Let h be the Fourier transform of (1 — P)y/. Then h(rf) G Ker(P) = 
Ker(N), and therefore h(g) is of the form 


c(rf)g + V <S> co(?]) + lo(?i) ®g, r) G R 4 . 


By taking the inverse Fourier transform, we see that (1 — P)y/ is of 
the form eg + d s u for some c G P'(R 4 ) and u G P'(R 4 ; A 1 ). Let y 2 G 
Cf > (M) satisfy y 2 (x) = 1 . Now (x, £) f WF(Py/) implies (x, £) f 
WF (/ — y 2 (eg + d s Lu)), since 


(37) 


x/ = PxS + (1 - -P)x/ = Pxf + eg + 


Suppose now that there are c G £'(M) and u> G £'(M ; A 1 ) such that 
(x, 0 £ WF(/ + cc/ + d s w). Then (x, 0 £ WF(P/) since 


yP(/ + eg + d s cn) = y Pf, 


and the pseudodifferential operator P does not move singularities. Fur¬ 
thermore, the formula ([37]) implies that (x, £) G WF(APy/) if and only 


if (x, £) G WF(Ay/). 


□ 
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10. Microlocal analysis of the light ray transform 


Let us now study the light ray transform in the limited angle case, 
that is, when U\ C M 3 is an arbitrary open set in Theorem |4.4[ Let 
y E Cq°(Mi). The restricted light ray transform yL belongs to the class 
of operators considered in [T4]. Here we will complement the results 
in [H] by computing the symbols of the restricted light ray transform 
and its normal operator with a suitable cutoff. 

It is well known that a restricted ray transform is a Fourier in¬ 
tegral operator associated with the canonical relation given by the 
twisted conormal bundle N* Z' of its point-line relation Z |15j . We 
will next give a parametrization of N*Z' and the symbol of yL in this 
parametrization. 

We take the Fourier transform with respect to y E M 3 and get 


TLf( V ,v) 


fim(s, y + sv)6 l 9 m e~ ir,y dsdy 
f lm (x)9 l 9 m e~ i ^ x '- x0v) dx. 


Here 9 = (l,u) and we have used the change of variables ( |23| ). We see 
that the Schwartz kernel of yL is the oscillatory integral 

(38) [ {x,y,v,rj)dr), x E M 4 , (y,v) 6l 3 x S 2 , 

Jr 3 


where a(x,y,v,rj) = (27r) 3 9 l 9 m and (ft(x,y,v,rj) = rj(y — x') + x°r]v. 
We write C = M 3 x S 2 . The critical set of (ft is 

C<t> = {(y,w,x-,-q) eCxI 4 x(l 3 \ 0 ); = 0 } 

= {(y,v,x,r})-, y = x’ - x°u}. 

The phase function (ft parametrizes the conormal bundle N*Z of the 
point-line relation 

Z = {(x' — x°v, v; x°, x') E C x M 4 ; v E S' 2 , x E M 4 } 

via the diffeomorphism C$ 3 (y,v,x,rj) i—X (y,v,<ft' ,x,(ft' x ) E N*Z. 
Note that the twisted conormal bundle 


N*Z' = {( y,v,rj,w ; a; 0 , x', £ 0 , O e ( T*C x T*R 4 ) \ 0; 

y = x - x°v, r) — w = x°g\ TvS a, £o = ~£v, 

(39) v E S’ 2 , r) el 3 , xE M 4 } 

is a canonical relation. Indeed, N*Z' C ( T*C \ 0) x (T*M 4 \ 0), since 
rj = 0 if and only if = 0, and — 0 implies £ 0 — —<£ 7 u = 0. 
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In the parametrization ( |38[ ) , the symbol of \L is 

(2^fl‘ x (y)e l 9”, 0 = (!,»), (y,v)€C. 


The symbol is of order zero (38) and hence yL is a Fourier integral 


operator of order —3/4 associated with the canonical relation N*Z 
see |18, Def. 3.2.2]. The theory of Fourier integral operators implies 
that \L gives a continuous map from £'(R 4 ) to T>'(C). Taking into 
account the conformal invariance, we have proven Lemma [4. 1| 

Let us next verify the claims in the end of Section [| about the sharp¬ 
ness of Theorem 14.41 


Lemma 10.1. Let (y,v,rj,w, x,£o,£') £ N*Z'. Then £ is lightlike if 
and only if C\t v s 2 — 0. Moreover, if f is lightlike then £' = —£ou. 

Proof. We have |£o| = |^| < |£'||r;| = |£ 7 | with the equality holding 
if and only if £' and v are collinear. But £' and v being collinear is 
equivalent with £,'\t v s 2 = 0. Thus £,'\t v s 2 = 0 if and only if |£ 0 | = |£ ; |- 
Suppose now that £ is lightlike. Then £ = cv for some c G M, and thus 

□ 


-c. 


£o = -?V = 

Lemma 10.2. Let (x, £) G T*M 4 \ 0. Then the relation C = N*Z' 
mapping T*M 4 to T*C satisfies the following: 

( 1 ) Iff is timelike, then C(x,f) = 0 . 

(2) If f is spacelike, then C(x,f) = {(x' — x°v,v,ff L x°f , \T v s 2 )'-, v G 
S'!}, where is the circle defined in Lemma\8.1 


that 


(3) Iff is lightlike, then C(x,f) = (x' + ^£', - J,£',0). 

Proof. If (■ y,v;r],w ) G T*C is in relation N*Z' with (x,f) G T* 
is, (■ y,v,r),w ) G N*Z'(x,f), then 

161 = \?v\ < |£|. 

Thus £ is lightlike or spacelike. Let us suppose that £ be spacclikc. We 
have y = x' — x°v, rj = f and w = x°f'\T v s 2 - Finally, £ 0 = — f'v is 
equivalent with v G Let us now suppose that £ in lightlike. Note 
that £o / 0 since otherwise £ = 0. Lemma 10.1 implies that v = — f'/fo 
and that w — 0. □ 

The restricted light ray transform yL maps singularites through the 
relation N*Z ', that is, WF ( x Lf) C N*Z' WF(/), see HH Th. 4.1.1, 
Th. 2.5.14], Hence Lemma 10.2 implies that WF(yL/) is empty if 
WF(/) contains only timclikc covectors. Let us now suppose that 
(x, £) G T*R 4 is spacelike and that the circle 

(40) = {x - x°v; v G S\} 
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does not intersect U\. Then Lemma 10.2 implies that WF(%L/) is 
empty if WF(/) is contained in a small conical neighborhood of of 
(x, £) in T*M 4 . 

The condition ( Jl2j ) is equivalent with S '}^ fl U\ ^ 0. Namely, if 
a null geodesic intersects x then, up to an affine reparametrization, 
it coincides with 7 (t]x,9) where 0 = (l,n) for some v G S' 2 . Now 
£(7(0; x, 6)) = £ 0 + €'v, whence v G 5} is equivalent with vanishing of 
£(7(0; x, 9)). Finally, 7 (t-,x } 9) intersects {0} x U\ for some v G S} if 
and only if 5'}^ fl U\ 7^ 0. 

Let us now suppose that (x, £) G T*M 4 
the null geodesic 

(41) h(t) = (x° + r, x' - t ^ _1 


\ 0 is lightlike, and consider 


0 o, 


r G 


Then N*Z' maps all the points (//(t),£) G T*M 4 , r G M, on the same 
point in T*C. We expect that this leads to artifacts when trying to 
recover lighlike singularities, however, we do not analyze this further 
in the present paper. 


11. Microlocal analysis of the normal operator 


Let us consider an operator A having kernel of the form (38), and 
use the following notation for the projections 


(42) 



When the projection p is an injective immersion, the diagram (42) is 
said to satisfy the Bolker condition. In this case the corresponding 
normal operator A*A is a pseudodifferential operator [T3j. In our case, 
the Bolker condition holds only outside the set 


Jzf = {(y,v,r),w; x,£) G T*C x T*M 4 ; £ is lightlike}. 


Lemma 11.1. The restriction of p on N*Z' \ A? is an injective im¬ 
mersion but dp fails to be injective on N* Z' fl . 


Proof. We may use (y,v,£',x°) gCx M 3+1 as coordinates on N*Z. In 
these coordinates 


P 


(y,v,£,x 0 f' | r „ S 2), 


dp = 


Ad 0 0 0 \ 

0 id 0 0 

0 0 id 0 

\0 * * ‘fks 2 / 
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where * denotes an element that does not play a role in the proof. We 


see that dp is injective if and only if C\t v s 2 0- By Lemma 10.1 this 
is equivalent with £ being lightlike. 

Let us now show global injectivity of the restriction. Let (y, v, £', ;r 0 ) 
and (y, v, £', x°) be in C x M 3+1 , and suppose that £ is not lightlike and 
that p maps the corresponding points in N*Z on the same point. Then 
y = y, v = v, £' = £', and x°£,'\t v s 2 = x 0 £'| t v s 2 - Moreover, x° = x° 
since £,'\t v s 2 0 as £ is not light like. □ 

We will show next that A*A is a pseudodifferential operator if we 
microlocalize away from the set Jzf. 

Lemma 11.2. Let (x,£) G T*R 4 \ 0. Then the composed relation 
C = ( N*Z ')~ 1 o N*Z' mapping T*M 4 to itself satisfies 

(1) //£ is timelike, then C(x,f) = 0. 

( 2 ) //£ is spacelike, then C(x,f) = (t,£). 

(3) If £ is lightlike, then C(x,f) = {(/r(r);£); r G M}, where p is 


the null geodesic (fl) 


Proof. The timelike case follows immediately from Lemma 10. 2| Let £ 
be spacelike, and let (y,w,7],w) G T*C be in relation N*Z' with (t,£) 
and with a point (x, £). We have £' = p = £'. This implies that 
also £° = — fv — £°. As £'|t„s 2 7 ^ 0 by Lemma 
x ( fi’\ Tv s 2 = w = t°£ / \t v s 2 implies x° = x°. Final 
x' — x°v, whence x! = x’. 

Let £ be lightlike. Then Lemma 


10.1 


y x 


the equation 
- x°v = y = 


10.2 


is the only point in T*C that is in re 
determines £ uniquely but x only up to a translation along p. 


says that (x 1 + f-£', —£', 0 ) 
ation N*Z' with (ar,£).° This 

□ 


Lemma 11.3. Let us define C = N*Z' \ TZ. The composition C 1 o C 
is clean in the sense of ca Th. 21.2.14], The projection H 2 from the 
intersection of (W 1 x C with T*M 4 x diag ( T*C ) x T*M 4 to (T*M 4 ) 2 is 
proper, and the fibers 7 r^ 1 ( p), p G (T*M 4 ) 2 , are connected. 

Proof. We write 

X = C~ l x C, Y = T*M 4 x diag (T*C) x T* M 4 . 

We need to show that X D Y is a smooth manifold and that 

(43) T P (X n Y) = T P X n T P Y , p G X n Y. 

Lemma |11.2 implies that X D Y coincides with the smooth manifold 
diag (C) after a reordering of the factors in the Cartesian product. We 
may use p = (y,v,f',x°) GCx M 3+1 as coordinates on X D Y. Notice 
that a tangent vector of A", parametrized by (5p, Sq) G T^ PjP fiC x M 3+1 ) 2 , 
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coincides with a tangent vector of Y if and only if dp dp = dpdq. This 


implies dp = dq since dp is injective on C. Thus (|43|) holds. 

H 4 ) 2 


A fiber 7r 2 1 ‘ 


(p), p G (T*M 4 ) 2 , is diffeomorphic with the circle S'/ if 
p = (x,£,x,£) and £ is spacelike, and it is empty otherwise. □ 

Lemma 11.3| says that if we cut away the lightlike covectors, then 
we can use the clean intersection calculus [ESI Th. 25.2.3], Let Xo e 
S 0 (T*M. 4 ) satisfy supp (xo) Ll = 0 where 

Jzfjr = {(#,£) ^ T*M 4 ; £ is lightlike}. 

Then A*Ax 0 is a Fourier integral operator with the canonical relation 
C~ l o C where C = N*Z'\S if, and Lemma 11.2 implies that C _1 oC is 
in A^(diag (M 4 ))'. Hence A*Ax 0 is a pseudodifferential operator. The 
calculus implies also that the order of A*Axo is —1 if the order of A 
is —3/4 since the fibers 7 t£ 1 (p) , p G (T*M 4 ) 2 , are one dimensional (i.e. 
the excess of the composition is one). 

If we are given the data Af we can not compute Axof, but we 
can choose pseudodifferential operators xo on M 4 and xi 011 C s ° that 
XiA = Xi^Xo modulo a smoothing operator. Altenatively, we could 
choose pseudodifferential operators xi and Xo, both on M 4 , so that 
XiA*Axo = X iA*A modulo a smoothing operator. By Lemma 11.2 
this follows if Xo = 1 in supp ( X i) and supp ( X o) hi = 0. However, 
we prefer using xi on C since we need also a cut off x 2 £ C'o °(Ui) due 
to the fact that we do not have data on M 3 \ li\. 


Proposition 11.4. Let x 1 € C^QR) satisfy supp (x 1 ) C (—1,1), and 
define 

Xi(v,v) = Xi(\v v \/\v\), V e S 2 , r) G M 3 . 

LetU\ C M 3 be open and let xi G C'^/Zfj). Then the principal symbol 
of the pseudodifferential operator L*xiX‘ 2 .L 


ob(z,£) 


(lei 2 - led 2 ) 1 / 2 



- x°v)9 j 9 k 9 l 9 m dv , 


where 9 = (1, v). 


Proof. We choose xo £ Cq°(— 1,1) such that xo = 1 in supp (xi), and 
define Xo(£) = Xo(|£o|/|£'|)- If (: y,v,r},w ; x,£ 0 ,£') G N*Z’ then rj = £' 
and £ 0 = -£'v = ~r)v. Thus \r]v\/\r]\ = |£o|/|£'| and xo(£) = 1 if 
Xi( v ,v) 7^ 0- Now Xi L = Xi-^Xo modulo a smoothing operator and 
supp (xo) hi 2^ = 0. Hence L*xiXiL is a pseudodifferential operator. 

The kernel of L*x.\XuL is of the form 

f e -i*(x,y,v£) e uK*,y™) a fa Vj rfidpdgdydv, x,ze M 4 , 

JCx I 3 xM 3 
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where a(z,v,rj) = (27r) 3 xi(u, v)X 2 {z' — z°v)9 j 9 k 9 l 9 m . We have 
</>(-, V, v, v) - <i>(x, y, v, $') = (■}] - C)y - 4>{z, v , r)) + 4>(x, v , £'), 

where 4>(x,v,t;') = ffx' — x°£,'v. Hence the kernel simplifies to 
5(rj - v, v )d v d£dv 


’ S 2 xR 3 xI 


D i£' (x'—z')—i (x l0 — z °) £' v 


' 5 2 xl 


a (z,v,g)dgdv 



0 i£.{x-z) * 


R 4 J S 2 


K*c5(v)a(z,v,£')dvd£ : 


where k^(v) = £ 0 + £'v. By (25) the symbol cr(z,£) of L*XiX 2 L i 


is 


a(^0 = (2vr) 4 (|n 2 -|eo| 2 )- 1/2 
2*Xi(M\?\) 


a (z,v,£')dv 


(lei 2 - leoi 2 ) 1 / 2 Jsi 


X 2 (z' - z°v)9 j 9 k 9 l 9 m dv. 


Thus a 0 (x,f) = cr(x , £). 


□ 


Proof of Theorem 4-4 ■ By the conformal invariance it is enough con¬ 
sider the Minkowski case. Let (x, £) be spacelike and choose the cutoff 
Xi in Proposition |11.4 so that Xi — 1 near |£o|/|£l- Then cr 0 (x, £) = 
c{tf)N ]klm , where c(£) ^ 0 and N^ klm is the tensor (28) with x( v ) — 
X 2 (x' — x°v). We recall that the visibility condition Jl2[ ) is equivalent 
to fl U\ 7 ^ 0, where the circle S ^ is defined by (140b . Thus we can 
choose non-negative y 2 G C™(Ui) so that x does not vanish identically 
on S^. 


Let P be the projection (35). Note that both N and P are homo¬ 
geneous of degree zero, and Wo = N + (1 — P) is elliptic near (x,£). 
indeed, if W 0 f = 0 then 

0 = PWof = NPf , 0 = (1 - P)W 0 f = (1 - P)f. 

Thus Pf G Ker(iV) = Ker(P) and / = Pf = P 2 f = 0. 

As Wo is the principal symbol of the zeroth order pseudodifferential 
operator 

W = C - 1 L* X iX2L + (l-P), 

there is a parametrix Q of W such that W — 1 modulo a smoothing 
operator near (x, £). Thus P = QWP = Qc~ 1 L*xiX 2 LP ■ Hence 

(x, 0 G WF (Pf) if and only if (x, £) G WF(L*XiX 2 PP/), 

and we apply Lemma [973 to conclude. □ 
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